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W)' Abstract 

' Tlie notoriously small X — a ^S"*" excitation energy of the BN diatomic has been calculated 

using high-order coupled cluster methods. Convergence has been established in both the 1-particle 
basis set and the coupled cluster expansion. Explicit inclusion of connected quadruple excitations 
T4 is required for even semiquantitative agreement with the limit value, while connected quintuple 
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excitations Ts still have an effect of about 60 cm ^. Still higher excitations only account for about 
10 cm^^. Inclusion of inner-shell correlation further reduces by about 60 cm~^ at the CCSDT, 



with earlier calculations and experiment, albeit with a smaller (and conservative) uncertainty. The 
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Q I and 85 cm^^ at the CCSDTQ level. Our best estimate, re=183±40 cm^^, is in excellent agreement 
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'O^i dissociation energy of BN(X ^U) is De=105.74±0.16 kcal/mol and Z)o=103.57±0.16 kcal/mol. 
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I. INTRODUCTION 



The lowest electronic excitation energy of the boron nitride diatomic is among the most 
vexing problems in small-molecule computational chemistry. Not only are the X and 
a ^S"'' states nearly degenerate, but the combination of moderate multireference character 
in the X state and pathological multireference character in the a state makes the 
transition energy Tl excessively sensitive to the electron correlation treatment. 

Martin et al.[l[, using multireference average coupled pair functional (ACPF) 
techniques!^, found the ^11 state to be the ground state and predicted Te=381±100 cm~^. 
These authors^ also found that (nowadays) commonly used coupled cluster methods such 
as CCSD(T) 4] yield qualitatively incorrect answers. In an MRD-CI study, predicted 
. Elaborate multireference calculations by Mawhinnev, Bruna, and Grein (MBG) 2], by 
PetersonPI, and by Bauschlicher and Partridge (BP)jy| obtained considerably lower Tg 
values of 241±160 cm'^, 190±100 and 180±110 cm'^, respectively. Wattsl?!, at the 
CCSDT (coupled cluster with all single, double, and triple excitations |^) level with a cc- 
pVQZ (correlation consistent polarized quadruple zeta[12|) basis set, found Te=844 cm~^, 
and conjectured that this serious overestimate was due to neglect of connected quadruple 
(T4) and higher excitations. Both Boese et al.P] and Tajti et al.jl^, in the context of 
high-accuracy computational chemistry protocols developed in their papers, found that, in 
strongly multireference systems, T4 can easily make energetic contributions on the order of 
the difference between the CCSDT and multireference values. (Denis^^ crudely estimated 
the effect of T4 by assuming error cancellation with higher-order T3 in the singlet but not 
the triplet state, and predicted Te=175 cm~^.) Finally, a very recent quantum Monte Carlo 
(QMC) study by Lu[l2] in the present journal found 178(83) cm~^, where the uncertainty 
band represents one standard deviation in the QMC approach. 

The two most reliable experimental estimates are the noble gas matrix IR measurements 
of Bondybey aud coworkersQ. 15-182 cm-, and the negative .on time-of-flight 
photoelectron spectroscopy value of Neumark and coworkers |l4i|. 158±36 cm ^. 

The purpose of the present work is to establish whether a converged result can be obtained 
at all from single-reference coupled cluster methods, whether this estimate is in agreement 
with the other theoretical approaches and experiment, and finally what is the breakdown of 
various contributions in the cluster expansion. 
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II. COMPUTATIONAL DETAILS 



All calculations were carried out using the general coupled cluster code MRCC of Kallay 
and coworkers [ifil. The Austin-Mainz version of ACES II 3] was used to generate the 
required integrals and molecular orbitals. Unless otherwise noted, the CCSDT/cc-pVQZ 



reference geometries of Watts [3] were used, re(X ^n) = 1.3302 and re(a ^S+) = 1.2769 A. 



n 



Correlation consistent |l7t] (cc-pVnZ), augmented correlation consistent |l8| (aug-cc- 
pVnZ), and core- valence" .elation co^.tentQ (ccpCV.Z, ba.. .e. were used 
throughout. The largest such basis sets used, cc-pV5Z, is of [6s5p4:d3f2glh] quality. Where 
appropriate, contributions were extrapolated to the 1-particle basis set limit using the 
A + B/L^ formula of Halkier et al. | " 



III. RESULTS AND DISCUSSION 



All computed values are given in Table I, compared with available experimental data. 

As expected, the CCSD results are grossly biased towards the triplet state (Te=4432 cm~^ 
at the basis set limit). Inclusion of T3 (connected triple excitations) is required for an even 
qualitatively correct result, although even the CCSDT basis set limit Te=827 cm~^ is 3-4 
times too large. Quasiperturbative T3 corrections such as CCSD(T) overcorrect, and wrongly 
predict a singlet ground state We conclude that CCSDT is the lowest acceptable level 
of theory for the reference geometry. Comparison of the CCSDT/cc-pVQZ and CCSDT/cc- 
pV5Z values suggests that the latter is converged to within 2-3 cm^^ with respect to the 
basis set. 

Inclusion of T4 (connected quadruple excitation) proved essential for anything 
approaching quantitative accuracy. The CCSDTQ/cc-pVQZ basis set calculations reported 
here involve 419 and 391 million amplitudes, respectively, for the singlet and triplet states. 
They ran for two weeks each on single AMD Opteron 846 processors. At the basis set limit, 
T4 reduces the transition energy by 514 cm~^. We do note — as we have previously notedj^ 
for other strongly multireference systems like C2{X ^S+) — that basis set convergence 
for the T4 contribution is fairly slow (unlike for systems dominated by a single reference 
determinant 9] ) . This can be rationalized in terms of very prominent double excitations in 
the singlet wavefunction: dynamical correlation relative to them will be dominated by double 
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excitations, which represent quadruple excitations relative to the reference determinant. At 
the CCSDTQ basis set limit, we obtain = 313 cm~^, in agreement with Ref.Q but still 
considerably higher than the other results. 

Connected quintuple excitations (T5) still reduce the excitation energy by about 70 cm~^. 
Comparison of the cc-pVDZ and cc-pVTZ results for this contribution suggests that it 
converges quite rapidly with the basis set. 

Connected sextuple excitations only affect Tg by -8 cm~^, while the contribution of still 
higher excitations was found to be negligible. Our best estimate for the valence-only FCI 
basis set limit is therefore Te=243±28 cm~^, where our error bar is the sum of all the 
amounts covered by extrapolations. 

Somewhat surprisingly, the effect of core-valence correlation is found to account for the 
discrepancy with earlier theoretical studies and experiment. At the CCSDT level, it lowers 
by 59±7 cm~^, while the differential T4 core- valence contribution reduces Tf, by an additional 
26±4 cm-i at the CCSDTQ level. (The all-electron CCSDTQ/cc-pCVTZ calculations, at 1 
billion amplitudes each, took about one day per iteration running OpenMP-parallel on four 
AMD Opteron 846 CPUs. Sub-microhartree convergence requires about twenty iterations. 
Our attempts to carry out CCSDT/cc-pCVQZ calculations met with failure for the triplet 
state. Because of the clearly erratic basis set convergence behavior of the CCSD(T) energy 
in this case, we have chosen not to use the larger basis set data at this level of theory.) 

Our final best estimate neglecting spin-orbit splitting thus becomes Te=158±40 cm~^, in 
excellent agreement with the earlier calculations (which likewise neglect spin-orbit splitting, 
it being almost an order of magnitude smaller than their stated uncertainties). Our error 
bar is probably somewhat conservative, as it assumes that no cancellation at all would occur 
between extrapolation errors in individual contributions. 

The spin-orbit coupling constant of the X state is calculated as -24.27 cm~^ at 
the CISD/cc-pVQZ (uncontracted, no o functions) level using MOLPRo|2l|. in excellent 
agreement with the experimental value|23| of -25.14 cm~^. Its inclusion pushes up both the 
present calculated value and all the earlier theoretical values by these amounts: our final 
best estimate thus becomes Te=183±40 cm~^. This agrees with the experimental value of 
Neumark and coworkers 11411 to within the respective uncertainties, and finds itself near the 

n 

upper edge of the interval given by Bondybey and coworkers |13l|. 

Finally, as a byproduct of this study, we obtain the dissociation energy of BN(X ^11) using 
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W4 theoryj23| as L>e=105.74±0.16 kcal/mol and L'o=103.57±0.16 kcal/mol (the uncertainty 
being a 95% confidence interval). This is somewhat higher than previous calculated D^. values 
of 105.2 kcal/mol^l] and 104.2 kcal/molQ. The zero-point vibrational energy (ZPVE) of 
2.17 kcal/mol was obtained by combining the accurate Uf. and UeXe for the singlet state 
from Ref.0 with the state difference in ZPVE from Ref.Q- In Ref.0, %TAE[(T)], the 
percentage of the total atomization energy resulting from (T), was proposed as an indicator 
for the importance of nondynamical correlation effects. We note that %TAE[(T)]=6.03% 
for the X ^Tl state (on the low end of moderate nondynamical correlation), compared to no 
less than 18.63% for the a state (among the most severe cases surveyed in Ref. j^]). 



IV. CONCLUSIONS 

Summing up, the notoriously small X ^11— a excitation energy of the BN diatomic has 
been calculated using high-order coupled cluster methods. Convergence has been established 
in both the 1-particle basis set and the coupled cluster expansion. Explicit inclusion of 
connected quadruple excitations T4 is required for even semiquantitative agreement with the 
limit value, while connected quintuple excitations T5 still have an effect of about 60 cm^^. 
Still higher excitations only account for about 10 cm~^. Inclusion of inner-shell correlation 
further reduces Tg by about 60 cm^^ at the CCSDT, and 85 cm~^ at the CCSDTQ level. 
Our best estimate, Te=183±40 cm~^, is in excellent agreement with earlier calculations and 
experiment, albeit with a smaller (and conservative) uncertainty. The dissociation energy 
of BN(X 3n) is De=105.74±0.16 kcal/mol and L'o=103.57±0.16 kcal/mol. 
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TABLE 1: X - a^J:+ transition ener gy (cm 
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valence correlation 
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827.4 


CCSDTQ-CCSDT 


-323.7 
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FCI-CCSDTQ56 


-0.9 












-0.9 


243.3 


inner shell corr. 




cc-pCVDZ 


cc-pCVTZ 


cc-pCVQZ 


cc-pCV5Z 






Best estimate Running total 


CCSD(T) 


-15.2 


-15.7 


-6.2 


+4.8 






+16.2 




CCSDT 


-36.0 


-52.4 










-59.3 


184.0 


CCSDTQ-CCSDT 


-12.1 


-21.6 










-25.6 


158.4 



Best estimate, this work 
Inch spin-orbif^ 
MRACPF, Martin et al.[l| 
MRDCI, Mawhinney et al.[2] 
ICMRCI, Peterson[5] 
ICMRCI, BP[6] 
QMC, Lu[12l 
Expt. (matrix) [l^J 
Expt.(gas phase) ^4] 



158+40"^ 
183+40 
381+100^* 
241+160^* 
190+100'* 
180+110^* 
178+83'* 
15-182 
158+36*= 



(a) cc-pVTZ basis set used on boron. 

(b) extrapolated from CCSD(T)/cc-pV5Z value and -189.1 cm"! at the CCSD(T)/cc-pV6Z level. 

(c) Expt. Ao=-25.14 cm-i ;2a|, calc. ^^=-24.3 cm'^ (this work). 

(d) Value does not include spin-orbit splitting in triplet state. 

(e) From observed To=0.031+0.004 eVll4t and ZPVE difference from Ref. [ij, assuming 4 cm ^ uncertainty on ZPVE 
difference. 
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